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Abslract 


Rcinfoicing earth with geosynthetics as a method of ground improvement has gained 
immense importance in leccnt times. The present study is concerned with two facets 
of icinforccd eaifch namely, foundation beds and retaining walls The foundation 
beds leinforced with multiple layers of geosynthetics are analysed using FEM The 
soil geotextile mtciface is modelled as a contact problem instead of the conventional 
techniciue of adopting spring elements Parametric studies have been conducted 
for a homogeneous linearly elastic granular soil to determine the effects of the 
number of layers of reinforcements, the natuie of load distribution and intensity, 
modulai ratio and Poisson’s ratio on settlements, tensile force distribution and 
stress distribution. The study reveals that the optimum placement depth is in 
the range of 0 4b - 0.5b and optimum spacing is 0 4b for all the cases. Lateral 
stress distribution indicates a strong effect of confinement while the shear stress 
distribution demonstrates the effectiveness of the reinforcement in reducing shear 
stresses in the soil mass The study has been validated by comparing results with 
FEM studies i elated to reinforced foundation beds. 

The study also deals with optimal cost design of geosynthetic reinforced earth 
retaining walls subjected to static and dynamic loading. Choosing the length and 
strength of the reinforcement as a feasible design vector the Sequential Uncon- 
strained Minimisation technique is used to arrive at an optimal cost of the reinforced 
earth wall. Optimal cost tables are presented for different combinations of loading 
and the developed procedure is validated by taking up an example problem. It is 
found from a typical example problem that savings of the order of 7 - 8% can be 
made over the standard design. 
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Chapter 1 


Introduction 

1.1 General 

Reinforcements with geosynthctics as a method of treatment of soils of lo^\ bear- 
ing capacity is now an established ground improvement technique. Although, the 
systematic study of reinforced earth did not start until recently, to the villager 
in tropical Africa and Southern Asia, the practise of building houses and roads 
on fibre-reinforced earth is an age-old art. Reinforced soil structuies are being 
used more frequently in system infrastructure such as bridge abutments, pavements, 
railway foundations, air strips, working platforms for oil drilling, etc. The concept 
of reinforced earth was developed in 1960’s by a French Engineer, Henri Vidal 
The concept involves the increase in load carrying capacity of a composite formed 
by the soil mass and reinfoi cement due to the development of tensile stresses in 
the reinforcement and the shear bond with the surrounding soil. Excellent field 
performance of the reinforcements confirm the advantages pertaining to this type of 
construction. 

Different tj-pes of Geosynthetics are available in the market, viz geotextiles, 
geogrids and geomats Action of each of these differs in many ways when used as 
an inclusion in the soil. In contrast to the geomats and geogrids, geotextiles do 
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noL offei an}'^ bending resistance Geognds offer more frictional resistance against 
pullout Geosynthetics are incieasingly being used in the construction of reinfoiced 
caith structures, such as, retaining walls, embankments, soil foundations, etc 

The construction of reinforced soil foundation has considerable potential as a 
cost-effective alternative to con\entional methods of support In this technique, one 
or more layers of geosynthetic reinforcement and controlled fill material are placed 
beneath the footing The purpose is to provide a suitable operating surface on which 
loads may be earned without the subgrade failing or deforming excessively 

Retaining walls as eaith structures are frequently constructed for a variety of 
applications, most common being bridge abutments and road construction They 
aic requued where a slope is uneconomical or technically infeasible When selecting 
a retaining wall type, mechanically stabilized earth (MSB) walls should always be 
considered MSB (i e. reinforced soil) walls are basically composed of some type of 
reinforcing elements, e. g gcosynthetics in the soil fill to resist lateral earth pressures 
When compared to conventional retaining wall systems, reinforced earth retaining 
walls have significant advantages They are very cost eff’ective, especially for walls in 
fill embankment cross sections Furthermore, these systems are more flexible than 
conventional earth retaining walls such as reinforced concrete-cantilever or gravity 
walls Therefore, they are suitable for sites with poor foundations and seismically 
active areas. 

A proper understanding of the performance of any reinforced structure is achieved 
through large scale model tests, However, economic considerations pertaining to 
the conduct of such tests have encouraged development of numerical and analytical 
modelling techniques Such techniques can be used to study, m detail, the various 
parameters associated with the mechanism and simulate varied field conditions. 

Sophisticated theoretical predictive models may be adopted for complementing 
conventional analysis on large/important projects where the anticipated conditions 



aic such that the validity of simplified approaches may be questioned 

In the following section, a brief literature re\iew pertaining to the behavior of 
a footing undei stiip loading resting on reinforced bed and the design aspects of 
leinforced carth-ietamirig walls is piesented Based on the literature review, scope 
of the woik has also been outlined, 


1.2 Brief Literature Review 

The behavior of footings resting on reinforced beds and reinforced earth walls have 
drawn the attention of geotechnical leseaich community and much woik has been 
reported in literature However, only those literature pertaining to strip loading 
on lemfoiced beds are presented here mainly due to the reasons of space and 
brevity Another reason for the decision to work only with this type of loading as 
major constructions like highways and embankments impose such type of loading 
analytical and experimental studies on reinforced earth walls has also been presented. 

1.2.1 Reinforced Foundation Beds 

In order to predict the behavior of the reinforced earth beds, the various methods 
that have been originally developed for unreinforced soils have further been extended 
to study the effect of inclusion The theoretical predictive techniques that have 
been pioposed, fall in one of the following four categories based on limit equilibrium 
methods, lumped parameter models, theory of elasticity solutions and finite element 
methods. A brief review of literature pertaining to each of the above methods is 
presented as follows. 

Limit Equilibrium Methods 

Giroud and Noiray (1981) made an analysis of the unpaved road behavior with 
and without geotextile as reinforcement. They considered the effect of geotextile in 
increasing the bearing capacity of the soil This method is based on the consideration 



of the tension membrane” effect for the case of plane strain reinforced unpaved road 
dcfoiming uiidei the action of a single application of dual wheel load A-ssuniing an 
allowable rut depth, load spread angle, suitable geometric configuration and firm 
anchorage of the reinforcement, the strain in the reinforcement was calculated. On 
knowing the strain, the tension in the reinforcement is computed Design charts 
weie proposed foi the reduction in road base thickness, due to reinforcement effect 
of geos 3 mtlietic, for various traffic intensities 

Ingold and Miller (1982) proposed a theory for bearing capacity of footings on 
icmforced clay considering an equivalent undramed strength of reinforced clay A 
number of layers of lemforcements spaced at a certain distance were considered 
Considering the force equihbiium of an element of soil spliced between two layers 
of reinforcement, they presented an expression for equivalent shear strength to the 
soil They derived an equation for the bearing capacity of reinforced cohesive soil. 
Comparison of the predictive results with the model test results has shown a close 
agreement. 

Lumped Parameter Models 

Madhav and Poorooshasb (1989) proposed anew model which consists of Paster- 
nak shear layers, Winkler springs and a newly proposed rough membrane to repre- 
sent the mechanical response of a granular fill geosynthetic-soft soil system. The 
displacements of soil and tension in fabric were obtained by satisfying the force 
equilibrium, coupled differential equations. They reported that the soil settlement 
decreases and tends to become uniform with increasing rigidity of shear layer It 
was proposed that the geosynthetic should be placed within the granular fill rather 
than laying it directly on soft ground. 

Poorooshasb (1991) developed a unified approach to the solution of three classes 
of problems, viz. (i) mats supported on point loads, uniformally distributed or 
symmetrically distributed loads, (ii) mats bridging over voids appearing m the 
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subgiade aftei coiistiuction. (in) mats placed oxer non-uniform forms of ground 
subsidence In the analysis, the reinforcements aie assumed to be rough enough so 
that no slippage takes place The granular fill is considered as a strain hardening 
plastic mateiial having a yield function, / The special assumption regarding mode 
of deformation is that all the vertical planes in unloaded system remain both vertical 
and plane after loading has been imposed The results demonstrated the usefulness 
of the transfoim function developed on the basis of the above assumptions, m 
obtaining solutions to the problems of geosynthetic reinforced granular mats 

Gofai and Bouidcau (1994) described an analytical model based on the stochastic 
theory of stress diffusion. The model accounts for the effect of fill compaction 
and confinement as well as the boundary interferences with the stress diffusion 
mechanism 

Ghosh and Madhav (1994) studied the settlement, confinement and membrane 
effects of reinforced granular fill- soft soil system These studies were undertaken 
by modifying the models proposed by Madhav and Poorooshasb (1988,1989) As in 
the previous study the Tough membrane’ element is incorporated for single layer of 
reinforcement assuming horizontal transfer at the soil-reinforcement interface. The 
model is generalized by incorporating non-lmear response of the soft soil and fill 
under plane strain loading condition. Parametric studies indicated that at large 
deformations, the vertical component of the tensile force in the reinforcement resists 
the applied load from acting directly over the soft soil. The membrane effect is 
significant for low values of shear stiffness of the granular fill, A reinforcement 
length of 3b xvas found to be sufficient to improve foundation response 

Ghosh and Madhav(l994) developed another model to incorporate the confine- 
ment effect of a single layer of reinforcement. It is quantified in terms of average 
increase in confining pressure due to the reinforcement, from which modified shear 
stiffnesses of the granular fill surrounding the reinforcement are obtained. The 
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parametric studies indicates that while the membrane action iinpro\es the footing 
i espouse, the confinement effect fuither enhances it The confinement effect is more 
pronounced when shear stiffness of the granular fill is large 

Shukla and Chandra (1994) proposed a foundation model to incorporate the 
compressibility of the granular fill by attaching a layer of Winkler springs to the 
Pasternak sheai layer The parametric studies indicates that the consideration of the 
compressibility of the granular fill results m a significant increase in the settlement 
of reinforced soil 

Theory of Elasticity Solutions 

Binquet and Lee (1975a) arc the pioneers in carrying out an anaKtical study on 
the bearing capacity of reinforced soil beds Using the Boussmesq’s stiess distri- 
bution, they calculated the shear stress distribution on the plane of reinforcement 
and the maximum tensile stress in the reinforcement They also defined the locus 
of the maximum shear stresses with depth The distance xq along the length of 
the reinfoicement where the shear stress is maximum, is a function of the depth 
of reinforcement and the soil properties They suggested values of maximum shear 
stress as 0 3go a-t a depth of 0.256 and around O.lgo at a depth of 26, where go is the 
intensity of the applied loading and b, the width of the footing. The normal force on 
the reinforcement due to the loading was obtained by integrating the vertical stress 
over the area of reinforcement. The authors considered three failure mechanisms 
viz. pullout, tension failure of the reinforcement and shear failure of the soil above 
the reinforcement. They assume that the tension m the reinforcement is inversely 
proportional to the number of reinforcing layers. Based on this assumption, they 
have given expression for the tension as a function of shear stresses and normal 
stresses at that depth, and the bearing capacity ratio (BCR), defined as the ratio of 
average contact pressure for the reinforced soil to that of unreinforced soil. 

Pitchuraani (1992) presented an analysis for the reduction in settlement due to an 
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inclusion of strip icmforcemcnt within the soil mass An elastic continuum approach 
\\as used and the solution given by Mindhn (1936) was numerically integrated to 
study tlie effect of shear stresses along the strip on the settlement reduction at the 
surface Curves presented show the effect of length of the strip, distance and depth 
of the strip on the settlement reduction. The optimum depth of placement of the 
topmost layer is found to be 0 256, b being the width of the footing 

Finite Element Method (FEM) 

Blown and Poulos (1981) demonstrated how a finite element model of reinforced 
eaith could be used to investigate the increase in beaiing capacity and stiffness of 
a foundation due to the placement of the reinfoi cement. The various components 
of the reinforced soil structure, are considered seperately. incorporating an elasto- 
plastic soil model obeying Mohr-Columb failure criterion The reinforcement strip 
is treated as elastic with zero flexural stiffness transmitting axial forces only. The 
model accounts for slip at the reinforcement-soil interface. Two cases, a footing on a 
homogeneous soil and a footing on a foundation incorporating a cavity or a pocket of 
soft soil were studied The results of the finite element study show similar patterns of 
behavior to the laboratory model test results. The analysis shows that the quantity 
of reinforcement required to produce some significant increase in bearing capacity 
is high, and since the limit state of the reinforcement soil bond is reached at an 
early stage, the reinforcement soil slip (rather than the stiffness of the reinforcing 
material) is the governing criterion. The analysis brings out that the provision of 
the reinfoicement helps in spreading the load, thereby causing mobilization of soil 
resistance over a wider area and a shallow depth. 

Andrawes et al , (1982) presented a finite element analysis to predict the load 
settlement behavior of soil geotextile system using variational approach to obtain 
the stiffness matrix for the soil elements. The geotextile was represented by linear 
elements with no bending stiffness, The soil geotextile interaction was assumed to 
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be only fiictional and simulated bj spniig elements of 2 ero length connecting the 
nodes of soil elements to those of geote\tile elements Results obtained fiom this 
anahsis were compared with those obtained fiom model tests which consisted of a 
geotextilc reinforced sand layer loaded to failure by a strip loading The measured 
and predicted data foi load settlement curve showed that, when geotextile was placed 
at a depth of 0 5 times the width of the footing, the difference in these values were 
within 10 percent of the peak footing load For the case where the geotextile was 
placed at a depth equal to the width of the footing, the difference is again within 10 
pel cent, but only upto 0 85 times the peak load Beyond this the predicted value 
dneiges and grossly overestimates the bearing capacity of the footing The reason 
could be the local shear failuie of soil, which could not be accommodated in finite 
element procedure. 

Love ct ab, (1987) formulated a finite element program, in order to verify their 
model tests capable of handling large displacements and strains The subgrade 
was modelled as elastic perfectly plastic material The fill material was modelled as 
an elastic frictional material obeying Matsuoka yield criterion. The reinforcement 
was treated as perfectly rough, to simulate no slip condition ^’lelding of the 
reinforcement was also not considered based on observation made in model tests. It 
was also assumed that the reinforcement cannot take compressive stresses. The 
results of this analysis compare well with those obtained with model tests for 
reinforced and unreinforced cases The results indicated that reinforcement prevents 
the shear stress on its upper surface from being transmitted to the clay layer below. 
The membrane action was mobilized at large displacements, leading to large increase 
in normal pressure both above and below the inclusion. The membrane tension 
profile displayed a slight dip directly underneath the footing, confirming that tension 
m reinforcement is not constant under the footing. 


9 



Rowe and Sodcrman(1987) discussed the choice of the finite element and constitu- 
ti\e models and validation of the results against bench mark solutions demonstrated, 
both fiorn field evidence and theoretical analysis, that the reinfoicement plays a 
lelativel}' small lole at low load levels since the soil is essentially elastic Significant 
strain m the geotcxtile begin to develop with increasing plasticity and in fact most 
of the strain develops aftei a contiguous plastic region has developed in the soil 
Burd and Hulsby (1989) presented a finite element model capable of computing 
the behavior of a leinforced unpaved road deforming m plane strain The calcula- 
tions wcie based on the use of Von hlises plasticity model for the clay subgrade 
and a model based on the Matsuoka yield function for the granular fill The 
reinforcement was modelled using elastic membrane elements The formulation was 
based on a large strain, large displacement approach. The mesh used was based 
on six-noded isoparametric triangular elements to model the soil and three noded 
membrane elements to model the reinforcement. The finite element method was 
used to back analyze laboratory scale model tests m which a footing has been jacked 
automatically into a model reinforced unpaved road There was a good agreement 
in the load displacement curve for reinforced and unremforced systems between the 
model test results and finite element results The authors commented that the finite 
element model may be used to obtain values for parameters that are necessary for 
a full understanding of the behavior of reinforced unpaved road, but are difficult 
to measure in physical tests. The authors presented plots of two such parameters 
VIZ strain in the reinforcement and normal stresses acting on the reinforcement. 
They concluded that the finite element model is useful in estimating values of those 
parameters that are difficult to measure through model tests 

Burd and Brocklehurst(1990) desnbed asmall displacement Finite Element study 
to investigate the mechanism of reinforcement in a plane strain reinforced unpaved 
road under the action of single monotonic load The fill and the clay are both 
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modelled by six-noded triangular elements with three noded membrane elements 
used to model the leinforcement which is placed directl} at the base of the fill over 
the whole w'ldtli of the mesh An elastic-perfectly plastic model was used to model 
the clay subgrade and an elastic-perfectly frictional model is used for the fill where 
the plastic behaviour is modelled by a non-associated flow rule The reinforcement 
is modelled as a perfectly elastic material with an infinite yield stress m tension and 
zero yield stress in compression. The only parameter varied was the reinforcement 
stiffness. No slip was allowed at the soil-reinforcement interface and geometric non- 
linearity was excluded The computed results agreed well with pieviously published 
limit state models A particular feature was that the increase in reinforcement 
stiffness causes substantial increase in the magnitude of the interface stress and 
reinforcement force. 

Patel and Chandreshekhran (1992) formulated a FEM program to analyze two 
structures, viz a strip footing and a road/pavement rail track on reinforced soil 
subgrades. The soils were represented by four noded quadrilateral elements and 
geotextile by straight linear elements. The interface was assumed to be purely 
frictional and simulated by spring elements of zero lengths connecting the nodes of 
soil elements to those of inclusion elements. The results indicated best improvement 
at a depth of 0.556, the depth of the top layer of the reinforcement. They also 
concluded that failure mass predicted by FEM analysis is deeper than that predicted 
by Bearing Capacity theory Varying F, and Eg values have shown that there is 
greater potential for improvement in weaker soils 

Wilson- Fahmy and Koetner (1993) presented a finite element model to represent 
soil-geotextile interaction using an incremental load transfer 1 ~ D finite element 
type of analysis to study the behavior of geogrids embedded in sand under a pullout 
condition. The conversion of 2 - D structure of geogrids into an equivalent 1- D 
structure while still retaining its 2 — D properties was elaborated 
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Ot(im et al , (1994) performed beaiing capacity analysis of geognd leinfoiced 
foundation ground using rigid plastic finite element method A senes of parametric 
studies by changing the reinforcement conditions, i.e , length, depth and strength 
of geognd were conducted Results showed that the bearing capacity increases as 
depth and length of reinfoi cement increases but there exists an optimum depth at 
which the maximum reinforcement effect is mobilized. 

Desai (1995) piesented a nonlinear finite element analysis of the problem of 
strip footing icsting on reinforced sand bed The non-linear behaviour of sand 
was modelled by a hypeibolic reversible stiess-strain relationship. The geotextile 
was modelled as a new isoparametric element It was discretized into a set of 
membrane and contact elements while the sand continuum was discretized into eight- 
noded isoparametric elements. They concluded that the geotextile reinforcement 
contributes to reduction in settlements of vertically loaded strips and the major 
principal stress direction at a point in sand is a function of magnitude 

Gharpure(1995) conducted FEM studies to predict the load-deformation be- 
haviour of a geotextile reinforced sand bed subjected to strip loading. The soil- 
geotextile interface was modelled by defining the contact conditions at the interface 
rather than the conventional joint elements The parametric studies were carried out 
to find the optimum placement depth of the reinforcement, variation of shear and 
tensile forces along the reinforcement, effects of the nature of load distribution, effect 
of modular ratio, etc. The author concluded that the optimum placement depth of 
a single layer of reinforcement is 0 6b and that at optimum placement depth, the 
modulus of elasticity of the geotextile if increased beyond 200 times the modulus of 
elasticity of soil, does not further reduce the settlements 

Raghvendra et ah, (1996) extended BLnquet and Lee’s approach for the design 
of reinforced soil bed as a two layered system. They used a general purpose Finite 
Element program, FEAP using four noded quadrilateral elements to find the optimal 
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thickness of the upper granular layer Stiess distribution obtained from FEM analy- 
sis was used to calculate non-dimensional parameters for estimating the mobilization 
of tension in the reinforcement They suggested that these non-dimensional paiam- 
eters together with selected design variables could be used to calculate the thickness 
and length of the reinforcement required to satisfy the performance criterion on 
settlement and bearing capacity 

Sithaiam et al,,(1996) conducted Finite Element Studies to obtain the nature of 
tensile distribution along the reinforcement for a reinforced soil bed A nonlinear 
confining stress-strain dependent behaviour model was used for the soil elements 
Finite Element Analysis was earned out with different types of reinforcement and 
loading was applied sequentially in increments. They concluded that the tensile 
force distribution along the length of the reinforcement is non-linear Higher the 
reinforcement stiffness, higher was the tensile force Another observation was that 
the location of the maximum tensile force depends on the type of the reinforcement. 

Experimental Investigations 

Bmquet and Lee (1975b) investigated the mechanisms and potential benefits 
of using reinforced earth slabs to improve the bearing capacity of granular soil. 
They conducted model tests with stiff footings on reinforced foundations for three 
conditions - homogeneous deep sand, sand above an extensive layer of very soft 
material and sand above a finite size pocket of very soft material. The reinforcing 
material consisted of 13mm wide strips of household aluminium foil. Tests were 
conducted for linear density ratio, defined as the total width of strips per unit width 
of footings, of 0.425. Three series of tests were conducted. The results of these tests 
showed that the load settlement and the ultimate bearing capacity of the footing 
could be improved by a factor of about two to four times above the bearing capacity 
of an unreinforced soil under identical conditions. Results were obtained regarding 
the effects of the depth of the top layer of reinforcement and number of layers of 
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rcinfol cement It was obser\ed that the bearing capacit}' ratio (BCR) continued 
to impiove with mcieasing number of lajers upto 6 to 8, beyond \\hich theie was 
little additional impiovement It i\as also observed that the optimum location of 
the top strip is at a depth of 0.31), b being the half width of footing and with further 
increase in this depth, the BCR value decreases. Thej’’ observed that pullout failure 
generally occurred with lightly reinforced slab {N < 2), whereas tie-breakmg which 
occurred m the uppermost layers was generally associated with heavily reinforced 
slabs (TV > 4) 

Akinmusiiru and Akinbolade (1981) conducted a senes of laboratory model tests 
with square footings on a deep homogeneous sand bed reinforced with flat strips of 
lope fibre material and undertook parametric studies. The results indicated that 
bearing capacity was maximum when the topmost layer of the reinforcement is at 
a depth of 0.5b from the base of the footing, where b is the width of the footing 
The optimum vertical spacing of the reinforcement was found to be less than 0.5b 
For spacing greater than that, the BCR decreased sharply and the rate of decrease 
became low for spacing greater than b. They concluded that increasing the number 
of layers beyond thiee does not contribute much to bearing capacity improvement 

Fragaszy and Lawton (1984) conducted a series of laboratory model studies 
to determine the influence of soil density and reinforcing strip length on the load 
settlement behavior of reinforced soil The tests were conducted on reinforced and 
unreinforced sand. The results showed that when BCR is calculated at a settlement 
equal to 10 percent of the footing width, the BCR is independent of soil density. 
When calculated at a settlement of 4 percent of the width of the footing, the 
percentage increase in the bearing capacity appeared to be less for loose sands 
than for dense sands. As strip length W'as increased from 3b to 7b, BCR increases 
rapidly but beyond that additional strip length does not appear to increase the BCR 
significantly. Comparison of their results with those reported hy Binquet and Lee 
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for eaith slabs sliowed good agreement in respect of the observation that the bottom 
layer of the leinforcemcnt breaks first The authors commented that the equations 
developed by Binquet and Lee are very sensitive to the value of friction angle at 
failure 

Guido ct al , (1986) presented a comparison of results of laboratory model 
tests to study the behaviour of geognd and geotextile reinforced earth slabs The 
parameters studied were the coefficient of friction between soil and geotextile, pullout 
resistance between geognd and soil, depth below the footing for the first layer of 
reinforcement, veitical spacing of the layers, number of layers, width size of square 
sheet of reinforcement and tensile strength of the reinforcement. The geotextile 
used was Du Pont Typar 3401 and geognd was Tensai SSI grid. They concluded 
that for both the cases, BCR decreased with increasing depth of top layer upto a 
critical depth of b, beyond which it remained practically constant As the vertical 
spacing was increased, BCR decreased. The number of reinforcement layers showed 
an optimum value of three, further increase in the number did not contribute 
significantly to BCR enhancement. 

Love et al , (1987) conducted small scale model tests to study the effectiveness 
of geogrid reinforcement, placed at the base of a layer of granular film on the surface 
of soft clay. To make the modelling as realistic as possible, the various components 
were scaled by factor of 4. The reinforcement used is a version, of Tensar SS geogrid. 
In the tests, monotonic loading was applied by a rigid footing under plane strain 
condition, to the surface of reinforced and unxeinfoxced system for a range of fill 
thicknesses and subgrade strengths Deformation of the subgrade and the geogrid 
were measured from photographs. Load in the footing was measured by a load 
cell and displacements by a Linear Variable Displacement Transducer (LVDT) 
The result of the model study indicated that the geogrid reinforcement tends to 
reduce the shear stress transmitted to the surface of the clay subgrade The authors 
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opined that the failure mechanisms m the clay are mobilized at small deformations 
of the fill and theiefore laige deformations are not necessary for the benefits of 
the reinforcements to be filled. At large deformations, where these are permissible, 
additional benefit was obtained from the membrane action of reinforcement 

Ramaswamy et al , (1992) conducted model tests on circular footings resting on 
reinforced clay they studied the influence of depth of a single layer below the footing, 
number of layers, thickness of grid material, anchorage length, sand-cover on both 
sides of the remfoicement and soaking of compacted soil on the BCR Similar to 
earlier studies, their conclusions were that bearing capacity increased with increase 
in the number of inclusion layers and thickness of geognd They concluded that 
there is substantial increase in BCR if the geognd reinforcement is covered by sand 
on both sides The bearing capacity of soaked reinforced soil was found to be greater 
than soaked unreinforced soil 

Giroud et al., (1993) with then laboratory shear box tests showed that the 
relationship between the interface shear strength and the stress normal to the 
interface is not linear (as it is assumed generally) and a hyperbolic expression 
represents a behavior close to the actual one 

Khing et al., (1993) conducted model tests on a strip foundation supported by 
a sand layer reinforced with layers of geogrids (punctured sheet drawn PPfHDPE 
co-polymer). The results indicated that the maximum benefits of geogrid reinforce- 
ments in increasing the bearing capacity was obtained when the depth of top layer 
of reinforcement was less than b, the foundation width. A reinforcement placed 
below tha foundation at a depth of more than 2,25 times the foundation width did 
not contribute to any increase in the BCR Maximum benefit was achieved when 
the minimum width of the geogrid layers was about six times the foundation width. 

Omar et al , (1993) earned out model tests for the ultimate bearing capacity of 
strip and square foundation supported by geognd reinforced sand They concluded 
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that foi de\clopnient of maximum bearing capacit)', the effective depth of the 
roinfoi cement should be about 25 for strip foundation and 1 45 for square foundation 
The maximum depth of placement of the first laj'er of reinforcement should be 
less than about b to take advantage of the intrusion The maximum width of the 
reinforcing layers required for mobilization of maximum BCR was found to be about 
85 for strip foundation and 4 55 foi square foundation 

Yctimogulu et al , (1994) performed laboratory model tests as well as Finite Ele- 
ment analysis to investigate the beating capacity of rectangular footing on sand bed 
reinforced with a single layer of reinforcement Both the methods revealed that there 
IS an optimum reinforcement embedment depth, 0 35, at which the bearing capacity 
was highest The studies also showed that there is an optimum vertical spacing 
of between 0 25 and 0 45 foi sand beds with multiple layers of reinforcement. The 
bearing capacity of reinforced sand was also found to increase with reinforcement 
layer number and size when the reinforcement is placed within a certain effective 
zone, which is approximately within 1.55 from both the base and edge of the footing. 

Adams et al., (1997) performed large-scale field tests to investigate the potential 
benefits of geosynthetic-reinforced soil Two different types of geosynthetics were 
used . stiff biaxial geognds and geocells. They found that BCR was maximum when 
the depth of the top layer was less than 0.55 At low strains, depth of the top layer 
within 0.255 showed good results Effective zone for the reinforcement was between 
0.255 and 0.55 for one and two layer case and between 0 5b and 1.55 for three layer 
case. 

1.2.2 Reinforced Earth Walls 

Investigations of reinforced earth retaining walls involve the state of stress inside 
the soil mass, the forces mobilized along the reinforcement, the forces acting on the 
wall facing and model tests on these walls Various design procedures have also been 
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developed for the reinforced earth walls The Literature review is presented in the 
following section 

Lee et ah, (1973) studied the use of reinforced earth for earth retaining structures 
The basic assumptions and results of the analysis were checked and found to be 
realistic by comparing the results with the observations made on model walls For 
ideal uncomplicated conditions, a cost analysis indicated that reinforced ear th walls 
may cost as little as one-half the cost of conventional walls Maximum tensile forces 
in the ties weie calculated and compared with the tensile strength to ensure adequate 
safety against breaking. Observations of failure heights for tie-breaking cases showed 
that there was little significant difference in failure heights for dense sand and loose 
sand They proposed that a better method of analysis of stress distribution within 
the ties and the backfill may lead to a more nearly optimum distribution of tie 
dimensions and therefore to more economical designs 

Kennedy et ah, (1980) presented an analysis to estimate the forces mobilized 
in the reinforced wall subjected to a surcharge strip load. An estimate of the 
maximum traction force in a reinforcing element was deduced by a modification of 
the Terzaghi’s formula for the conventional rigid retaining wall They also proposed a 
simple distribution for the traction force in a reinforcing element, and to estimate the 
force transmitted to the facing of the wall. They also concluded that the surcharge 
strip loads applied outside the top extremity of the failure wedge do not generate 
significant forces either in the reinforcement or on the facing 

Jewell and Milligan (1985) considered deformations of soil structures reinforced 
by extensible materials to ensure that the deflections at the end of construction and 
due to creep of the Reinforcement during the life of the structures are of a serviceable 
magnitude.They presented an analysis procedure linking an equilibrium stress field 
in which a constant mobilised angle of friction is assumed, with a strain increment 
or displacement field in which a constant angle of dilation is assumed Charts were 
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derived foi horizontal deflection at face of the wall and the vertical settlement behind 
the crest of the wall for an ideal length lavout at ideal and uniform spacing Similar 
charts were also derived for truncated length layout 

Labba and Kennedy (1986) conducted both experimental and theoretical studies 
to assess the maximum tensile forces mobilized in a reinforced wall subject to 
a vertical surcharge strip load or the combined action of vertical and horizontal 
sui charge strip loads They developed a simple design method for determining the 
maximum magnitude of the tensile force and its distribution with depth of the 
leinforced eaith backfill The ability of the wall system to transfer stiess from 
overstressed legions to those legions where the reinforcing elements not yet reached 
their full fiictional or strength capacity, in order to maintain internal equilibrium, 
was incoipoiated in the analysis. It was concluded that the stress transfer was 
greatly influenced by the surcharge load magnitude and its location, i e distance 
from the wall face 

Juran et ah, (1989) conducted laboratory model study on the performance 
and failure mechanisms of reinforced soil retaining walls, using different reinforcing 
materials, namely woven polyester geotextile strips, plastic grids, and non-woven 
geotextile strips. The model walls were instrumented to obtain measurement of 
stresses in the reinforcements, displacements at different points along the reinforce- 
ments, displacement of the facing, and to identify, using coloured sand, the failure 
surface in the soil. The results indicated the importance of confinement of the 
reinforcement on the structural performance, 

Leshchinsky and Boedeker(1989) presented an approach for stability analysis of 
geosynthetic reinforced earth over firm foundations. The internal stability analysis 
was based on variational limiting equilibrium and satisfied all equilibrium require- 
ments Two extreme inclinations of reinforcement tensile resistance were investi- 
gated: orthogonal to the radius defining the geosynthetic sheet, and horizontal, 
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signifying the as-installed position It was found that the required embedment length 
IS longei for the orthogonal inclinations The external stability analysis was carried 
out by using the bilinear wedge method The results of both type of analysis were 
piesentcd in the form of design charts 

Juran et al ,(1990) argued that neither of the working stress method and the limit 
equilibrium method consider the fundamental requirement of strain compatibility 
between the soil and the reinforcement nor do they provide any evaluation of the 
effect of the soil dilatancy and extensibility of the reinforcement on the mobilized 
tension foices and structural stability. They presented a strain compatibility design 
approach that is fundamentally based on the analogy between the plane strain shear 
mechanism that develops along a potential failure surface in the actual structure and 
the response of the reinforced soil material to simple shearing. 

Claybourn and Wu (1993) reviewed six published design methods, namely. Forest 
Service method, Broms method, Collin method, Bonaparte method, Leshchmsky 
and Perry method, and Schemertman method They made a comparative study 
of the design concepts involved in each of these methods and the results obtained 
by each of these methods. The computed results were then compared with the 
observations made on the test walls. They found that the differences in the design 
results stem primarily from the significant disparity in defining allowable reinforced 
strength and safety factors, and to a lesser degree, from discrepancies among the 
methods of analysis viz the lateral earth pressure distribution. 

Singh and Basudhar (1993) presented the application of a generalized approach 
to the estimation of the lower bound bearing capacity of reinforced soil retaining 
walls by using Finite Element technique in conjunction with non-linear programming 
to isolate optimal solution. The analysis was based on a rigid-plastic model for re- 
inforced soil, treating it as a macroscopically homogeneous anisotropic material No 
slipping phenomenon was considered at the soil-reinforcement interface The stress 
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field was tlie design vector subject to constraints arising from interface equilibrium 
and external boundary conditions as well as no-yield criterion The results when 
compared with data available in Literature showed good agreement 

McGown et al , (1994) outlined the Limit State approach to the design of 
polymei reinforced walls, slopes and embankments. The approach was shown to 
be applicable and it was suggested that a single integrated approach based on 
identifiable kinematic models can be generated for the reinforced structures 

Wong et al , (1994) conducted a senes of model tests to study the failure model 
of a geotextile reinforced soil wall The backfill was uniform fine sand and the 
reinforcement was a woven polyester with a strength of AlKnfvn A uniform 
surcharge pressure of up to 250 KPa was applied on top of the model A two part 
wedge type of failuie with the slope plane passing between two reinforced layers 
was checked in. the design. They concluded that sliding failure may be critical if 
surcharge load is applied behind the reinforced block. Their results did not support 
pullout and overturning failure modes. 

Karpurapu et al., (1995) described Finite Element models that were used to 
simulate the behaviour of two instrumented large-scale geosynthetic reinforced model 
walls The walls were constructed using dense sand- fill and layers of geosynthetics 
attached to two different facing treatments. The model walls were taken to collapse 
using a series of uniform surcharge load applied at the sand-fill surface. The authors 
adopted a modified form of constitutive hyperbolic model that included a dilation 
parameter to simulate the behaviour of granular soil The results showed that FEM 
models and constitutive models can predict important features of wall performance 
fairly accurately 

Bastick and Segrestin (1996) applied double wedge equilibrium method to rein- 
forced earth and analysed taking safety coefficients or weighing factors into account. 
They compared this method with the commonly used coherent gravity method and 
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g.abil stability method Consistent results were shown pro\ided that the safet\ uas 
correctly ensuicd 

Hjodo et al , (1996) studied the displacement of a reinforced earth wall using 
a small model involving two kinds of geosynthetic strips namely rubber and high 
density polyethylene (HDPE) Varying lengths of strip leinforcement were tested in 
dry Aio Sand. They observed that higher the tensile stiffness of the reinforcement, 
smaller is the wall displacement for full transfer of active force to the reinforcement 
Furthermore, the strains measured along the strip provided information on the 
behaviour of the strip when subjected to, firstly active force and later an additional 
pull out force 

Poibaha and Goodings (1996) tested 24 reduced scale models of vertical and 
steeply sloping (1H:6V) reinforced soil walls using Kaolin as backfill, reinforced with 
a non woven geotextile simulant and loaded to failure under increasing self weight 
in the geotechnical centrifuge. The models were constructed on either firm or rigid 
foundation and different lengths of reinforcements were used. The tests showed 
that in vertical walls, failure developed entirely within the reinforced zone nhen 
LjH > 0.75. For steeply sloped walls, this threshold occurred when LfH > 0.67 
Better performance was observed for models on firm foundations Tension cracks 
in the backfill led to stress concentration in geosynthetics and it was suggested to 
avoid their development in practice. 

Ling et ah, (1996) conducted a seismic analysis of reinforced soil structures using 
a pseudo-static approach. Based on the analysis, design charts were developed for 
a reinforced wall and the effects of seismic inertia force were investigated. For soils 
with small friction angle, or for large seismic intensities, direct sliding stability 
concerns lead to impractically long reinforcement lengths. Therefore, the authors 
proposed that a permanent displacement limit should be used in the design so that 
direct sliding from an earthquake can be tolerable. The authors obtained yield 
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accleration of the reinfoiced soil block in a two-part wedge mechanism for a factor 
of safety against direct sliding equal to unity 

Nishimara ct ah, (1996) investigated the resistance against earthquake of ten 
geogrid-rcinfoiced soil walls following the Hyogo earthquake They found that 
although the earthquake caused some settlements and cracks in the foundation, 
the walls themselves were free of deformations Analysis based on two design 
methods proposed by GRB and PVVRI in Japan, revealed that stability was critical 
at horizontal seismic coefficients of 0 1 to 0.2 The failure modes were circular slip or 
two-part wedge sliding They proposed that stability will be enhanced by adequately 
increasing the length of the gcogrid in the upper portions of the reinforced zone 

1.3 Motivation and Scope of the work 

Considering the importance of the use of geosynthetics as a reinforcing material 
for foundation beds, much work has been done by researchers all over the world. 
The various approaches discussed above for the analysis of reinforced foundation 
beds have their own pros and cons The Limit Equilibrium approach, albeit simple, 
suffers from the limitation of not considering the deformations of the soil. Similarly, 
in theory of elasticity approach, the analysis approach is simplified by considering 
an equivalent elastic modulus for the reinforced zone. It is also not able to predict 
the nature and magnitude of stresses at the critical interface of soil and geosyn- 
thetic. The Lumped Parameter models are more general and can be applied to a 
variety of situations. In terms of computational cost and understanding the overall 
behaviour of the system, these models have an advantage But their utility is again 
limited if stresses and displacements at various points in the continuum are to be 
found out Model tests, though quite reliable, have economic constraints since it 
becomes extremely difficult and expensive to perform these tests mth all different 
combination of parameters. 
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FEM modelling enjoys a considerable advantage over all other approaches and 
experimental observations, since parameters can be easily varied for a compielien- 
sive study and details of stress distributions and displacements can be determined 
throughout the system This is particularly useful for studying the soil-geosynthetic 
interface. 

Review of literature reveals that when multiple layers of reinforcements are 
used, the reinforced zone is generally homogenized and analysis is earned out b\ 
consideiing the zone to be verj' stiff, having very high modulus of elasticity It has 
already been pointed out that such procedures do not give the correct representation 
of the stresses developed in the medium. As such theie is a need to analyze the 
reinforced bed considering each reinforcement as a sepeiate entity. Gharpure(1995) 
developed a computer program utilizing FEM to analyze reinforced sand beds with 
a single layer of reinforcement. The objective of the present work is to improve 
upon the developed computer program in order to incorporate multiple layers of 
reinforcements. 

The importance of reinforced earth retaining walls is easily recognizable. Al- 
though a number of design approaches have been developed and analytical studies 
earned out, optimal cost studies to check the economic feasibility of such walls 
have not yet been attempted .As such an attempt has been made here to develop a 
general procedure to arrive at an optimal cost estimate of a retaining wall reinforced 
with geosynthetics formulating the problem in the frame of non-linear programming. 
Studies have been made with regard to optimum cost of both geogrid and geotextile- 
wrap walls under static and dynamic loading. 

Chapter 2 deals with the Finite Element Analysis of the reinforced foundation 
beds. The optimum cost design of the geosynthetic reinforced earth walls is discussed 
in Chapter 3 In Chapter 4, the scope of the future work has been presented. 
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Chapter 2 

Finite Element Analysis of 
Reinforced beds 

2.1 Problem Formulation 

2.1.1 General 

Finite Element Method (FEM) is well recognized as a very powerful and versatile 
numerical tool used for solving continuum mechanics problems Numerous examples 
of its application to model reinforced soil can be found in the literature. As already 
discussed in chapter 1, the same has been adopted m the thesis also. In the following 
sections, the soil-geotextile interaction mechanism is described and then the use of 
the FEM in finding the response of reinforced soil bed subjected to superimposed 
loads has been presented. The procedure was originally developed by Gharpure 
(1995) for a single layer of reinforcement. The same has been extended here to 
study the effects of multiple layers of reinforcements. For the sake of completeness, 
the analysis procedure is presented in detail 
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2.1.2 Soil- Geotextile Interaction Mechanism 

Remfoiced Foundation beds are classified as internally stabilized load bearing sys- 
tems which aie dependent upon the interaction mechanism of the soil and rein- 
forcement Within this system, shear transfer occurs through the interface by 
the mobilization of the induced tensile force m the embedded reinforcement The 
process by which forces and strains are transmitted from the soil to the reinforcement 
depends on the soil-remforcement friction in case of granular soils, soil-remforcement 
interlock in case of geogrids and adhesion m case of undrained loading of soft clays. 
The tensile force distribution within the leinforcement solely depends on interface 
stiess transfer law ( Mohi-Coulumb’s Law). The tensile forces induce corresponding 
compression in the soil as long as there is no slippage between the soil and the 
reinforcement, This additional compression increases the confining pressure on 
the soil. As a result, the soil can take more axial load. This is the ‘confinement 
effect’ Since mobilized shear stresses confine the surrounding soil, the effect will be 
maximum directly below the footing. Refer Fig 2,1. 

Moreover, at large strains, the lateral displacements of the granular fill due 
to the applied load are large enough to mobilize the full shear resistance at soil- 
reinforcement interface. As a consequence, the geotextile, which does not have any 
bending strength, takes the shape of a hanging string. Due to the tension induced, 
it lifts the overlying soil. The applied load is distributed over a wider area on the 
soil underneath Vertical pressure and settlement m the zone below the footing are 
thus reduced Additionally, differential settlements are reduced due to induction of 
extra pressure and settlement on the soil away from the footing. Refer Fig. 2.2. 
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2.1.3 Statement of The Problem 

Fig 2 13 shows a strip load of width b having (a) uniform intensity of loading (for 
flexible foundations) and (6) paiabolic distribution of loading (for rigid foundations) 
acting on a semi infinite elastic half space A number of reinforcing geosynthetic 
ia\ers are placed from a depth d with spacing s. The settlement response of the 
footing and the stress behaviour is to be studied for different reinforcement and soil 
conditions The soil is assumed to be homogeneous with constant Young’s Modulus. 

The method developed has been validated by taking up case studies and compar- 
ing the obtained results with experimental and theoretical solutions available in the 
hteiature. Paiametric studies have been conducted to study the effect of impoitant 
factors on the stiess-deformation behavior of reinforced beds. 

t 

2.2 Analysis 

2.2.1 Assumptions : 

The following assumptions are made in the analysis : 

1 The soil is considered as a linearly elastic material 

2 The geo textile is also considered as a linearly elastic material 

3. Length of the footing is very long in comparison to the width of footing so 
that plane strain condition prevails. 

4. The Soil-Geotextile interaction is based on 

• Friction between soil and geotextile, 

• Adhesion between soil and geotextile, 

• String effect. 
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5 In compai isoii with the soil dimensions, the thickness of geotextile is negligible 
However, thickness is considered while calculating stresses in geotextile 


2.2.2 Finite Element Formulation of The Soil Bed : 


The finite element procedure used is basically that foi solving continuum load- 
displacement problem. It consists of four mam steps discretization of the contin- 
uum, selection of approximating function; derivation of elemental stiffness matiices 
and analysis of the element assemblage. 

The Problem domain is discretized by dividing it into a number of elements. The 
soil has been discretized using four noded rectangular elements and the geotextile 
inclusion has been discretized using linear elements. The complete problem geometry 
and mesh details are shown m Fig 2 3 

For a typical soil element, the following approximation to the dependent variables 
u and V (displacements in x and y direction respectively) is chosen 
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which can be written as 

M = (2.2) 

Here the functions (shape functions) are bilinear functions of natural coordinates 
((, 77 ) and {u}® is the elemental displacement vector. The chosen shape functions 
satisfy the completeness and compatibility conditions required for convergence. For 
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a plane strain pioblem, the strain vector has only three components 


w = 


The strain vector is related to the displacement vector {it} by 

[e] = W{w} 

where 

1-^ 0 1 

dx ^ 

|il = 0 i 

J- A. 

dy dx _ 

Substitution of approximation (Eqn. 2 2) into (Eqn 2 4) leads to 

{«} = (£1[W1{«}‘ = lB|{u)» 

where the matrix (B) is given by 
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For a linearly elastic material, the stress-strain relation is usually expressed as 

[^1 = 1 D ]{£} (2 8) 

where the elasticity matrix [D] for a plane strain problem is given by 


[D] = 


(1 + i/s)(l “ 2Us) 
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0 0 (1 - 2v,)/2 
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Tlic potential encigy H of an elastic body can be expressed as 


n = i - IjufindA - liu^Mdl (2.10) 

where {«}, {e}, {a}, {/} and {t} are vectors of displacements, strains, stresses, 
body forces and surface tractions respectively denotes the area of the body and 
r is the boundary of the body where the tractions are specified, Substitution of the 
relations (2.2), (2 6), (2 8) into the expression (2.10) for the potential energy leads 
to 

n = E - {x}"'!?}’! (2-11) 

e=l ■* 

I-Iere n is the number of elements, 

[kf = [ \Bf\D][B]dA (2 12) 

IS the element stiffness matrix and 

{<!}' = / [Nfif)dA + i (Wf (2.13) 

JQc 

is the element load vector. Further, is the area of the element e and is that 
part of its boundary which is common to F , 

The expression (2.11) can be wiittcn in terms of the global matrices as follows : 

n = Imnmu] - {(/^{Q} (2.14) 

and the minimization results in 

[JiTJM = {Q}. (2 15) 

Here 

IK] = t W' 

e=l 
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IS the global stiffness matrix, 


{Q} = (2 17) 

IS the global load vector and {U} is the global displacement vector The summation 
sign in equations (2,16), (2.17) means the matrices and {qY are to be expanded 
first to the global size and then added. 


Note that the evaluation of B-matnx (Eqn. 2 7) involves x and y derivatives of the 
shape functions A^,. However, A/, are given as the functions of the natural coordinates 
The derivatives can be calculated as follows By chain lule 
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The matiix [J] is called the Jacobian matrix, It can be evaluated from the relation- 
ship between the two coordinate systems The relationship is given by 

* = i:W.(C.>?W (2-20) 

t-l 


v = 'tN,(c,n)yt ( 2 . 21 ) 

1=1 

where (iCj,yj) are the nodal coordinates and Ni are the same shape functions as 
used in the approximation of u and v. Thus the chosen element is an isoparametric 
element. 
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2.2.3 Finite Element Formulation Of the Geotextile : 


The Geotcxtile layer is repiesented by straight line elements without any bending 
stiffness as shown in Fig 2 3. For the two different types of mechanisms, the 
geotcxtile element is modeled differently 

For the confining effect of geotcxtile due to tension, the element is modeled as 
an axial element with linear approximation for the displacement in x direction The 
same vaiiational approach as used for 2-D soil elements is used The corresponding 
stiffness matrix becomes 


■ 



-1 
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( 2 . 22 ) 


where Ag is the cross sectional area and is the length of the element. 


For the string effect, a transversely distributed load p as shown m the Fig. 2.4 is 
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Figure 2A: Definition Sketch for String Effect 


considered. The variational functional for this case is 
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wlicic V IS the displacement in v~dircction, T is the tension in the stimg and S is 
the vertical component of T at each end Again using the same procedure as used 
for 2-D soil element, the minimization of vaiiational function gives the following 
stiffness matrix 



1 -1 

-1 1 


(2 24) 


Note that the stiffness matrix for the confinement effect is related to the x- 
direction displacements while that due to the string effect is related to the y-direction 
displacements The combined stiffness matrix can be arranged as 


^9 -^9 

0 

~^ 9^9 

0 

0 

r 

0 

-T 

~^g^9 

0 

Eg-Ag 

0 

0 

-T 

0 

T 


(2.25) 


Finite Element Modeling of the Interface : 


In general practice, the soil geotextile interface behavior is modeled by introducing 
interface elements. This can be achieved in a number of ways, like, the use of 
joint element of zero or nonzero thickness, node compatibility spring element or by 
composite layer model 

The joint clement formulation is derived on the basis of relative nodal displace- 
ments of the soil elements surrounding the interface as shown in Fig. 2.5, where the 
subsciipts n and s indicate the normal and tangential (shear) components while the 
subscript r stands for the relative quantity. 

Another common approach is to join the adjoining three nodes by springs as shown 
in Fig, 2.6 In this type of formulation, the spring elements ensure compatibility of 
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c = normal siress at interface 
n 

X = shear stress at interface 
k = normal stiffness of interface 

n 

k = shear stiffness of interface 
s 

= relative normal displacement 
= relative shear displacement 


Figure 2 5‘ Joint Element of Zero Thickness 

the displacements between the adjoining nodes until the Mohr- Coulomb’s criteria is 
reached. It replaces the compatibility conditions by failure conditions and dilatancy 
equation once the interface strength is exceeded based on some assumed constitutive 
relationship {e.g Andrawes et al (1983) used a hyperbolic model) to represent the 
interface behavior 

The composite layer model involves a beam-merabrane-type four noded element 
with large deformation capability. This provides flexural stiffness and local mem- 
brane effect but does not allow slip between the reinforcement and soil. Therefore 
this is particularly suitable for geogrid type of reinforcement. 

In the present work, no interface element is used but the interface is modeled as a 
contact problem. Referring to Fig. 2.7, the soil-geotextile interaction is represented 
by the following contact conditions 

1. Force equilibrium conditions : For vertical forces at contact nodes 

(n).+ W). + (n)l. = 0 (2 26) 


’ = K 


Interface Thickness=0 
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and for horizontal forces at contact nodes 


{F,)r + {F,)j+{F,)k = Q. (2 27) 


2. Displacement conditions No separation is allowed between the soil and 
geotextile hence the vertical displacement condition at contact nodes is defined 
as 

(u), = (v)j = (v)jt ( 2 . 28 ) 

Horizontal displacement condition depends on the occurrence of slip between 
the soil and geotextile at a particular node Foi no slip, horizontal displace- 
ment compatibility is defined separately for the top and bottom soil nodes as 
slip can occur on either side of reinforcement. Thus 


(«), = (u)^. 

(2.29) 

(^i)j = {u)k 

(2 30) 


For slip between the soil and geotextile, nodal force condition describing the 
Molir-Coulomb criterion replaces the displacement compatibility condition as 
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Figure 2 7 Modeling of the interface 



follows. 

(^x), = {Fy%tan4>* + cT, (2.31) 

{Fa:)j = iFy)jtanf + c*l^ (2 32) 

where <f>* is the frictional angle between the soil and geot extile and c* is 
adhesion per unit area For parametric study adhesion is taken as zero. 


2.2.4 Analysis Of Element Assemblage : 

After evaluating the stiffness matrices for the individual elements, the stiffness 
matrix [/-f] for the whole system is assembled. The basis of assembly is that the 
value of the displacement at a node is the same for each element sharing the node. 
As the reinforcing layer nodes are not connected with the soil nodes, the assemblage 
takes the shape as shown in Fig 2.8 


The next step is to modify the stiffness matrix, (a) for Soil-Geotextile intci action 
mechanism as explained in the previous section and (b) for boundary conditions 
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Figure 2 8. Assemblage of Stiffness Matrices 








which includes zero oi nonzero displacements at specific boundary nodes (The roller 
supports show that only vertical displacements are allowed and fixed support shows 
no displacement m eithei direction is allowed). Then the load vector {Q} is calcu- 
lated from the given loading condition (uniform, parabolic etc.). This is followed by 
the solution of the simultaneous equations for the unknown displacements. Fmally 
from the known displacements, the nodal forces and strains and stresses for each 
element aie calculated to yield the complete solution. 

To chcclt the slip conditions at interface, the forces at these nodes should be 
known Hence the solution requires iterations In the first iteration, no slip is 
considered between the interface nodes and the check for slip condition is made in 
the second iteiation The iterations continue till slip conditions for all the nodes for 
two successive iterations remain the same 
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2.2 Results, Discussion and Conclusions 

2.2.1 Input Parameters 

The footing of width b is placed on a granular bed . The size of the domain for FEM 
analysis is chosen to be 5.5 times the width of the footing from the central axis along 
X direction and 3 times the width of the footing along Z direction Owing to the 
symmetry of the problem, only half of the layer need be anahzed The geosynthetic 
layers are placed such that these extend by 1 5 times the width of the footing beyond 
the edge of the strip loading. Analysis have been cained out with meshes of size 
24 X 16 elements and 12 X 8 elements. 

Soil Parameteis . 

Angle of Internal Fiiction, 4> = 35°; 

Poisson’s Ratio, i/ = 0 49; 

Modulus of Elasticity of Soil, = 30 MPa; 

Geosynthetic characteristics : 

Modulus of Elasticity of Geotextile, Eg — 4 GPa, 
thickness = 7mm; 

2.2.2 Parametric studies 

Four cases have been investigated. 

1. Unreinforced case This has been simulated using identical modulus of elas- 
ticity values for soil and geotextile 

2 Single Layer of reinforcement. 

3. Two layers of reinforcement. 

4. Three layers of reinforcement 
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Table 2 1 Effect of number of lajers on settlements Ijniform Distribution of loading 


Case 

Maximum 

Settlement 

(cm) 

percent 

reduction 

unremforced 

4 42 

0 0 

Single layer 

4.07 

7 87 

two layer 

3 86 

12 61 

three layer 

3 70 

16 23 


The results obtained for these cases with homogeneous soil subjected to uniform 
distribution of loading are as follows: 

I Settlernent profile: 

The settlement profile for the four cases are shown m Fig. 2 9 As expected, at 
z = 0, settlements along X direction are least when three lajers of reinforcements are 
placed Settlement is observed to reduce more and more with addition of each layer 
of reinforcement .To get a quantitative idea, the corresponding percentage reduction 
in maximum settlement over the unreinforced case is illustrated in Table 2 1 The 
Fig. 2 9 also shows that increase in the number of layers of reinforcements causes 
significant reduction in settlements upto a distance of 0.56 from the centre of footing 
along X direction Beyond this distance, there is no appreciable impact of increasing 
the number of reinforcing layers on the settlement characteristics 

To get an idea of the effect of nature of loading on the settlement profile, the 
settlement profile was plotted for the case of parabolic distribution of loading also. 
The profile obtained as in Fig. 2.10 shows similar trends to the uniform loading 
case. The percentage reduction in maximum settlement over the unreinforced case 
is illustrated in Table 2.2. Expectedly, the maximum settlements for the Parabolic 
loading case are more than the corresponding ones for the Uniform loading case. For 
a comparative view, the settlement profile for the uoreinforced case and three layers 
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Figure 2.11: Comparison of surface Settlement Profile for UDL and PDL case(unreinforced and 
three layer case) 


Table 2 2’ Effect of number of layers on settlements’ Parabolic Distribution of loading 


Case 

Maximum 

Settlement 

(cm) 

percent 

reduction 

un reinforced 

4 86 S 

0.0 

single layer 

4 48 

7.72 

two layer 

4 26 

12.31 

three layer 

4 09 

15 83 
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of reinforcements are plotted for both the loading distributions in Fig, 2 11. The 
fig shows that a rigid footing, coriesponding to the parabolic loading case, settles 
mote at the centre of the loading as compared to the flexible footing Be\ond the 
edge of the footing, XfB = 0 5, the parabolic loading case shows lesser settlements 
than the Unifoim loading case. 

I Effect of embedment on settlements: 

To study the effect of embedment depth ( the depth at which the first reinforce- 
ment layer is placed ), the embedment ratio, d/b has been varied from 0 3 to 0 8 with 
increments of 0 1 for all three cases of reinforcements The obtained results have been 
presented in Fig 2.12 It can be seen from the figure that as the embedment ratio is 
varied from 0.3 to 0.5, settlement reduces, but any further increase in the embedment 
ratio causes increased settlements. Table 2 3 show's maximum settlements 

for vaiious d/b ratios and the percentage reduction in settlements for the corre- 
sponding value obtained for unreinforced case It can be noted that the percentage 
reductions in settlement over the unremforced case for different d/b ratios fall within 
a very narrow zone for a particular case The optimum placement depth for all 
reinforcement cases is 0.5b. It agrees excellently with experimental investigations 
earned out by Akinmusuru et al(1981) and Ramaswamy et al(1992) The maximum 
percentage reductions in settlement values at optimum placement depth are 7.87, 
12.61 and 16.3 for single layer, two layer and three layer case respectively 

Similar studies have been conducted for a Parabolic distribution of loading. The 
results obtained for the three cases are presented in Fig. 2 13 For all the cases, 
the maximum settlements show a decreasing trend as d/b is varied from 0 3 to 
0 4. Further increase of embedment ratio causes an increase of settlements too 
Table 2,4 shows maximum settlements for various d/b ratios and the percentage 
reduction in settlements for the corresponding value obtained for the unreinforced 
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Table 2 3. Effect of embedment ratio on settlements: UDL case 


case 

d/6 

maximum settlement 
(cm) 

percent 

reduction 

unreinforced 

00 

4 4205 

00 


0.3 

4 1159 

6 89 


04 

4 0806 

7.69 

single 

0.5 

4 0727 

7.87 

layer 

0.6 

4 0809 

7.68 


0.7 

4 0975 

7.31 


08 

4.1181 

6.84 


0.3 

3 9072 

11.61 


04 

3.8676 

12.51 

two 

05 

3 8631 

12 61 

layer 

06 

3 8788 

12 25 


07 

3.9050 

11.66 


08 

3.9362 

10.96 


0.3 

3 737 

15.54 


04 

3 6999 

16 30 

three 

05 

3 7030 

16.23 

layer 

06 

3.7270 

15.69 


07 

3.7617 

14.90 


08 

3 8009 

14 02 











case. Thus, the optimum placement depth for Parabolic Distribution of Loading for 
all cases of reinforcements is 0 45, although the variations in percentage reduction 
m settlements are almost insignificant for a particular case of reinforcement. The 
maximum percentage reductions m settlement over the unreinforced case at the 
optimum embedment ratio are 7.72, 12.3 and 15 84 for single layer, two layers and 
three layers of geosynthetic respectively. 

I Effect of spacing on settlement: 

To study the effect of spacing between the layeis at the optimum placement depth, 
the spacing ratio, s/6, has been varied from 0.2 to 1.2 for two layer case and 0.2 to 
0 8 for three layer case It is evident from Fig 2.14, that the settlements reduce as 
s/b is varied from 0 2 to 0 4 and then show an increasing trend as s/b is increased 
further. The variation is more pronounced for three layer case. This agrees well 
with model tests conducted by Yetimogulu et al(1994) and Singh (1988) 

Quantitatively, Table 2 5 shows maximum settlements for different s/b ratios and 
percentage reductions over the unreinforced case It cleaily shows that the optimum 
spacing for both two and three layer case is 0.4b and the coi responding percentage 
reduction m settlements being 11 57 and 16 27 respectively, 

B Effect of Poisson’s Ratio on settlements: 

The effect of Poisson s ratio on settlements at the centre of the loading has been 
studied and the results are shown in Fig 2.15, The Fig 2 15 shows that an increase 
in Poisson’s ratio over the whole range of varying from 0 3 to 0 49 (diained and 
undrained condition) leads to an increase m percentage reduction in settlement over 
the unremforced case. 

The effect is more pronounced for the undrained case (i/, = 0 49) as is evident 
from Table 3 6 The corresponding maximum percentage reduction in settlement 
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Table 2 4 Effect of embedment ratio on saclements PDL case 


case 

d/b 

maximum settlement 
('em) 

percent 

reduction 

unremfoiced 

00 

4 837 

00 


0 3 

4 51 

7 14 


0.4 

4 482 

7 72 

single 

05 

4 489 

7 57 

layer 

06 


7 21 


0 7 

4 531 

6.71 


08 

4 551 

63 


0 3 


12 04 


04 

4 259 

12 30 

two 

05 

4 275 

11 98 

layer 

06 


11 36 


0 7 

4 34 

10 66 


0 8 

4 377 

9 88 


03 


15 64 


04 


15.84 

three 

05 

4.114 

15.3 

layer 

06 

4 154 

14.47 


07 

4.199 

13 55 


08 

4 24 

12 72 


Table 2 5 Effect of spacing ratio on settlements, {djb 0 5) 


case 

s/6 

maximum settlement 
(cm) 

percent 

reduction 

unreinforced 

EQ 

4.4205 

00 


on 


11.48 




11.57 

two 



1130 

layer 


3 9409 

10.85 


1.0 

3 9643 

10 32 


1.2 

3.9894 

9 75 


0.2 

3.7149 

15 96 

three 

0 4 


16 27 

layer 

IH 

3 7445 

15 29 


IIB 


1 13 82 
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Table 2 6 Effect of Poisson’s ratio on Setliemcnts, (d/b) = 05 



unremforced 

one layei 

tuo layer 

three layer 



settl 

settl 

% 

settl 



% 



(cm) 

(cm) 

diff 

(cm) 



diff 


1 

919 

1 766 

-7 97 

1 657 


1 573 

-18 04 

0 45 

2 

615 

2 481 

-5.12 

2 383 

-8 89 


-11.92 


3 

128 


-3 67 

2 958 

-5 44 

2 896 

-7.42 

0 35 

3 

574 


-1 83 

3 454 

-3 34 


-4 64 

0 30 

3 

986 

3 943 

-1 09 

3 903 

-2 08 

3 867 

-2 99 


values aie 7 97, 13 66 and 18 04 for single layer, two layer and three layer case 
respectively The reason for such behaviour is that m the undrained case, the soil 
IS at its stiffest state and so the setllcments are less than those for drained case 


I Effect of Modular Ratio on settlements 

The Modular Ratio is the ratio of the stiffness of the geotextile and the stiffness of 
the soil and thus, reflects the strength of the reinforcement. Fig 2 16 shows the effect 
of Modular Ratio, Eg/Es, at d/b = 0,5, on the maximum settlements at the centre 
of the loading for the three reinforced cases The figure shows that as the Modular 
Ratio increases, i. e. stiffer reinforcement is chosen, the settlements decrease at an 
increasing rate, but the decrease is at a constant rate beyond Modular Ratio equal 
to 400. 


I Variation of tensile force in reinforcement: 

Due to the soil-reinforcement interaction, tensile forces are mobilized in the re- 
inforcement. The tensile force in the top layer of reinforcement have been plotted 
in Fig. 2,17 for all the three reinforced cases. Predictably, the tensile forces are 
maximum for the single layer case and minimum for the three layer case There 
is a sharp decrease in tensile force from the centre to the edge of the loaded area 
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Thereafter, the tension in the geosjnthetic vanes moie gradually till it becomes zero 
at the edge of the reinfoi cement 

For the two layer case, upto a distance of 0 36 from the centre in the X direction, 
the top layer of the geosynthetic shows more tension than the lower layer (Fig 2 18) 
Beyond that distance, the lower iajer shows greater tension The results are compa- 
rable qualitatively to the FEM studies conducted by Burd and Brocklehurst (1990) 
and Love et al (1987) Model tests conducted by Love (1987) and Gourc et al 
(1985) also show a dip in tension profile directly beneath the footing The extent of 
significant tension is an indication of the required width of the reinforcement in the 
soil bed 

Tensile force distribution has also been plotted for the case of Parabolic Distri- 
bution of loading in Fig 2 19 As in the previous case, the mobilized tensile forces 
are maximum for the single layer case A comparative plot of tensile forces in the 
first layer of reinforcement for the cases of Uniform and Parabolic case of loading 
shown in Fig 2 20 predictably shows a direct relation between the tensile force and 
the load intensity. The increase in tensile force is in the range of 10-15% for the 
parabolic loading case over the uniform loading case, directly below the footing. 
Beyond the edge of the footing, the Uniform loading case shows greater tension in 
the reinforcement for both single layer and three layers case 

An attempt was made to study the effect of load intensities and Modular Ratio 
on the tensile force distribution along the reinforcement Fig 2 21 shows the tensile 
force distribution for load intensities of 10, 100, 200 and 500 KNfm for a Modular 
ratio of 200. It illustrates that the reinforcement becomes more effective at higher 
load levels as greater mobilization of tensile forces takes place at these levels. 

To study the effect of stiffness of the reinforcement on the tensile force distri- 
bution, tensile forces are plotted in Fig. 2.22 for a constant loading of IQQKN/m 
and varying modular ratio It can be observed that a stiffer reinforcement shows a 
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distinct peak tensile stress neai the centre of the loading ^^hlle the reinforcements 
with ^ery low stiffness show a uniform tensile force distribution The results of the 
stud}' of effect of load intensities and the modular ratio on tensile force distribution 
agree well qualitatively with similar FEM studies conducted by Sitharam et al 
(1996) 
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I Distribution of stresses: 

Normal stress: 

Plots of distribution of vertical normal stress at different vertical cross-sections 
aie taken at 0 063m, 0.19m, 0 31m, 0 44m, 0 94m, 1.6m, 2.81m and 5 11m from the 
centre of the footing Figures 2 23 to 2 26 show the distribution of vertical normal 
stress for all the four cases Expectedly, foi sections nearer the centre of the loading, 
stresses are higher than for those further away from the centre The stresses follow 
the general trend of Theory of Elasticity solution Stresses tend to decrease sharply 
for the sections nearei the centre wheieas the variation is more gradual for sections 
further away Similar trends are observed for all the four cases 

In Figures 2 27 and 2 28 the stresses for all the cases at each vertical section 
are compared For shallow depths, the stresses are nearly the same for each case 
but for depths greater than the reinforcement depth, the three layer case shows 
more stresses than the two layer, single layer and unreinforced case The percentage 
increase over the unreinforced case varies m the range of 7-10 for three layer case, 
6-8 for two layer case and 2-4 for single layer case For sections at 1.06m and 1.6m 
which are at distances greater than b from the centre of loading, the three layer case 
shows the least stresses for shallow depths The percentage reduction being of the 
order of 3-4 % 

The tension developed in the reinforcement due to the mterfacial friction induces 
a horizontal force in the opposite direction in the soil mass which acts as a lateral 
restraint, producing the effect of confinement. The lateral restraint produces an 
additional vertical force on the effective area of the inclusion which causes additional 
normal stress in the reinforced zone. This may be the cause of increased vertical 
stresses in the top layers as observed from Fig. 2 27 and Fig. 2 28. Perhaps, due 
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to reduction in the confinement effect from tlie centre, sections at 1 06m and 
beyond show lesser stresses for the reinforced cases 

In the available literature, such studies regarding stresses using the confinement 
effect have not been presented and thus, the values obtained by the present study 
could not be compared with the other solutions The above observation is contrary 
to the general expectation that increasing the number of layers of reinforcements 
causes reduction in the stresses 

Raghvendra et al (1996) studied the effect of increased number of reinforcing 
layers on stresses by considering the reinforcement zone being stiffer for more number 
of reinforcement layers. They have chosen the stiffness of the reinforced zone as ten 
times that of the zone below it The results obtained by the present study, after 
modifying the developed computer program on similar lines, agree qualitatively with 
those obtained by Raghvendra et al from FEM studies (Fig 2.29 and Fig. 2.30) 
Due to inadequacy of data available for the problem solved by the authors, only 
qualitative comparison could be made Since all other results of parametric study 
are consistent with experimental and theontical studies in this field, the reason for 
the apparent anomalous behaviour may he m the ‘confinement effect’ Even though 
care has been taken to get convergent solutions with appropirate mesh sizes, the 
calculation of stresses may not be that accurate as for that one needs a very fine 
mesh Effect of finer mesh on the stresses needs to be investigated further. 

For the two layer case, a plot of normal stress in the reinforced zone agrees, qual- 
itatively, with the Finite Element studies conducted by Burd and Houlsby (1989). 
The plot, Fig 2.31 indicates reduction in normal stress across the reinforcement. 
The stress difference arises due to the curvature of the reinforcement and tends to 
reduce the magnitude of the normal stress acting on the soil below the geotextile. 

Later cd stress: 

The benefits obtained from reinforced soil results from the generation of frictional 
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forces at soil-reinforcement interface These forces are analogous to the increase in 
confining pressure, thus restricting lateral strains in the soil As a result, Fig 2 34 
and Fig. 2 35 indicate higher lateral stresses for the reinforced cases as compared 
to the unreinforced case For sections beyond the reinforced zone m X direction, 
the diminishing of confinement leads to higher lateral stresses for the unreinforced 
case Figures 2 32 and 2 33 show the distribution of lateral stresses, individually, 
for the unreinforced case and the three layer case. At distances beyond the edge of 
the footing (A = 1 6m, 2 81m, etc), where the tensile forces reduce, the reduction 
in confinement effect causes a corresponding reduction in lateral stresses for the 
reinforced cases As a result, at these sections, the lateral stresses are less for the 
reinforced cases as compared to the unreinforced case 

Shear Stress: 

The lateral streses in the fill generate shear stresses. If the reinforcement is ab- 
sent, these shear stresses are sustained directly by the soil and so have a detrimental 
effect on the bearing capacity In the reinforced case, the shear stresses are assumed 
to be sustained wholly by the geosynthetic Thus, the presence of the reinforcement 
causes a reduction in the outward acting shear stresses leading to better performance 
of the fill under the superimposed load. Plots of shear stress for the four cases are 
shown in Figures 2 36 to 2 39 The comparison, as in previous cases, is shown 
in Fig. 2 40 and Fig 2.41. The plots highlight the effectiveness of increasing the 
number of reinforcing layers in reducing the shear stresses. The fact that the shesr 
stresses are sustained wholly by the reinforcement is brought out by observing the 
abrupt reduction in shear stresses in the reinforced zone for all the three cases 
Shear stresses acting on the reinforced zone are plotted for the two layer case. 
The resulte agree qualitatively with the findings of Burd and Brocklehurst (1990). 
They have showed (Fig 2 42) that the shear stresses reached a peak at about 0 76 
from the centre of the loading The present study shows (Fig 2 43) that the shear 
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stresses reached a peak at a distance of about 0 56, the edge of the loaded area, from 
the centre of the footing 
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Figure 2 23 Vertical Distribution of Normal Stresses Unreinforced case 



Figure 2v24 Vertical Distnhution of Normal Stresses* single layer case 
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^ Figure 2.3Q: Vertical Distribution of Normal Stresses (Validation problem, BrlE, = 10) 
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Figure 2 34. Vertical Distribution of Lateral Stresses Comparison 
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Figure 2 39 Vertical Distribution of Shear Stresses Three layer case 
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Figure 2 41 Vertical Distribution of Shear Stresses- Comparison 
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2.2.3 Conclusions 


The general conclusions that can be drawn from the present study are presented as 
follows 

• The computer program originally developed for a single layer of reinforcement, 
when extended to multiple layers of reinforcements, functioned well and the 
predicted results are found to be compatible with experimental and theoretical 
studies 

• Increase in the number of layers of reinforcements causes reduction in settle- 
ments, the reduction being of the order of 15 - 16% for a three layer case as 
compared to unreinforced case 

• The optimum placement depth for uniformly distributed strip loading is 0 56 
for all the cases of reinforcements For Parabolic distribution of loading, the 
optimum placement depth is 0.56 for a single layer case and 0 46 for two and 
three layer case 

• The optimum spacing of the reinforcements is 0.46 for both two and three layer 
case 

• Settlement reduction for all cases is most effective when soil is in undrained 
state 

• As Modular ratio, .Bg/E'a increases, settlements reduce at an mcreasing rate, 
but for EgjEi > 400, the decrease is at an almost constant rate The trend is 
similar for one, two and three layers of reinforcements 

• Tensile forces in the top layer of the reinforcement decrease with increase in 
the number of reinforcement layers Increase in the strip load on the footing 
causes greater mobilization of the tensile forces in the reinforcement 
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• Shear stresses are reduced sharply in the reinforced zone, the reduction being 
greatest for the three layer case 

• Due to confinement effect, lateral stresses show an increase Avith increase in 
the number of layers of the reinforcement 
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Chapter 3 


Optimal Cost Analysis of 
Reinforced Earth Walls 


3.1 Introduction 

The use of Reinforced earth walls for a variety of purposes was highlighted in 
Chapter 1 . Their importance is well recognized. The use of geogrids or geotextiles 
rather than metallic strips (ties) is a further development of the Reinforced earth 
concept. Geosynthetics offer a viable and often very economical alternative to 
metallic reinforcements for both permanent and temporary walls, especially under 
certain environmental conditions The Chapter deals with an optimum cost design 
of Geosynthetic reinforced earth walls. 

3.1.1 Brief description of the wall 

A sketch of the wall is shown in Fig. 3.1. The main item is the soil backfill which 
is placed as the wall is constructed. For the problem considered herein, the backfill 
IS limited to cohesionless, free draining material, i.e. sand , having the following 
properties : density 7 and angle of internal friction if>. The backfill is strengthened 
by reinforcing members, which in the present case are geogrids and geotextiles. Since 
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system for long-term applications. The facings, additionally, enhance the aesthetic 
appeal of the completed wall For geogrid wall. Modular Concrete units (MCU) are 
the most common facing used They are popular because of their aesthetic appeal, 
easy availability and relatively low cost Geotextiles are more commonly used m 
wrap- around-faced MSE walls Such facings are used for temporary structures, 
for walls that are subject to significant post-construction settlements and where 
aesthetic requirements are low 

3.1.2 Design approaches 

A number of approaches have been developed for the design of geosynthetic soil 
walls. Among them, the Limit Equilibrium method and the Working Stress de- 
sign method are more common The Working Stress method basically rely upon 
restrictive assumptions with legard to the state of stress in the soil. The Limit 
Equilibrium method essentially use conventional slope stability analysis, modified 
to account for the reinforcement effect, for the global stability of the reinforced 
soil mass For retaining walls, all the methods normally considered applicable 
to routine design, use limiting equilibrium analysis to determine factors of safety 
against failure, however, they vary in their assumptions regarding stress distribution, 
failuie surfaces, safety factors and the inclination of the reinforcement at the failure 
surface The various design methods used in practice are Forest Service method 
(1977), Broms inethod(1978), CoUm method(1986), Bonoparte et al. mettiod(1987), 
Leshchmsky and Perry method(1987), Schmertmann et al method(1990) These 
methods typically presume a planar failure surface through the reinforced mass 
described by the classical Rankine active failure condition. The reinforcements 
extend beyond the assumed failure surface and are considered to be tension resistant 
tiebacks for the assumed failure wedge. As a result, these methods are collectively 
referred to as tied-back wedge analysis methods 
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For calculating maxiinum stress developed m the ties or reinforcements, three 
approaches are possible These are designated as Rankine method, Coulomb 
Force method and Coulomb Moment method. The Rankine method considers the 
equihbiium of a single representative critical element of soil at any depth, whereas 
the Coulomb methods consider the overall stability of the entire wall. All these 
methods predict maximum tension in the lowermost layer of the reinforcement But, 
the Rankine method necessitates a minimum embedment depth beyond the failure 
plane for all ties or geosynthetics, while the Coulomb method does not require that 
all the ties should extend behind the failure plane Hence, the Coulomb method 
computes the shortest length of the reinforcement , whereas the Rankine method 
leads to the longest tie. Thus, the Rankine method is more conservative 

The earliest method to determine the lateral earth pressure on a retaining struc- 
ture subjected to dynamic loading was developed by Mononobe(1929) and Ok- 
abe(1926). The approach is basically a pseudo-static approach in which the Coulomb’s 
Earth Pressure equation was modified to account for the additional vertical and 
horizontal forces induced during an earthquake. The M-0 method was developed 
for a dry, cohesionless backfill. 

Factor of safety is defined as the ratio of the forces or moments resisting instabil- 
ity to the forces or moments tending to cause instability. Full mobilization of the soil 
strength is assumed to occur. This approach is unlike the one used for evaluating 
slope stability wherein the factor of safety is used to reduce strength parameters. 

Most methods available for designing reinforced soil walls utilize the same general 
input parameters . First, the wall geometry must be defined in terms of height and 
face inclination. Most methods allow for a surcharge behind the top of the wall face. 
Secondly, the strength properties of the reinforced soil mass must be described The 
parameters used typically consist of soil unit weight and internal angle of friction, 
reinforcement strength and parameters describing the frictional interaction between 
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the soil and the reinforcement 

Review of the literature revealed that eventhough lot of studies have been under- 
taken to develop methodologies for the analysis of reinforced earth walls, no work 
has been undertaken to develop methods for their optimum design Such an attempt 
has been made here 

3.1.3 Statement of the Problem 

Fig 3 2 shows the cross-section of a reinforced earth retaining wall. Given the 
wall geometry, the soil properties and the reinforcement properties, the problem 
is to define the length and stiength of the reinforcement that gives the least cost 
estimate for construction. 

3.2 Analysis 

Optimal cost design of retaining walls when reinforced with either geogrids or 
geotextiles is earned out Standard design procedures developed by Federal Highway 
Authority (FHWA) have been adopted This particular analysis procedure is based 
on the Rankine Approach The same is cast in the frame of optimal design procedure 
in the following sections. 

3.2.1 Assumptions 

• Wall should be perfectly flexible and frictionless 

• When walls are inclined, it should not interfere in the formation of the outer 
rupture plane and the soil should not slide on the wall 

• As per Rankine’s theory, the lateral pressure against the wall increases linearly 
with depth. 
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. The failure wedge slopes at an angle (e = 4b + ^1) to the horizontal. 


3.2.2 Design guidelines 
I Input parameters 

The input parameters for the optimal design procedure are presented in Table 3 1. 
Due to the lack of available manufacturer’s data, the values of the allowable tensile 
strength of the geosynthetic are interpolated 

I Design variables 

The design vectors chosen for the formuJation are related to the reinforcement 
They are -0(1) = /=Length of the Reinforcement D(2) = Ta=Long-term allowable 
strength of the reinforcement. 


I Design Steps 

The Active Earth Pressure coefficients for both internal and external stability were 
computed using the Rankine Earth Pressure Theory. For backfill at an angle /?, the 
coefficient of active earth pressure for external stability is defined as, 


= cos X 


cos — v'cos^ ^ — cos^ I 


cos /? + V'CQS''^ i3 - COS^ (pk 
The coefficient of active earth pressure for internal stability is defined as, 


(3.1) 


— COS fix 


COS P — ycos^^— COS^^ 


(3 2) 


cos (3 + P — cos^ <l)f 

The minimum recommended wall embedment depth is 0 45m Therefore, the 
design height of the wall is incremented by the addition of the minimum embedment 
depth of the wall. 
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Table 3 1 Input Parameters for Design of MSB wall 


Parameter 

Symbol 

Value 

Pleight of Wall 

h 

3 — lOm 

Minimum Embedment depth of wall 

^dmin 

0 45rri 

Angle of Internal friction of the fill 

<i>f 

35° 

Unit wt of the fill 

If 

20KN/m? 

Backfill inclination 


0.0 

Anglo of Internal Friction of the backfill 


30° 

Unit wt. of the backfill 

lb 

ISKNlnF 

Allowable Tensile Strength of the 

T 

** Q 

30 ~ 6(}KN/m 

Geosynthetic 



Minimum spacing between the 


0.2771 

reinforcements 



Minimum Embedment length of the 


1.0m 

reinforcement 



Design Factors of Safety: 



Overturning 


2.0 

Sliding 

FSst, 

1.5 

Beaiing Capacity 

FStc, 

2.0 

Reinfoi cement Strength 


1.5 

Reinforcement Pullout 

FSpij 

2.0 

Cost Factors’ 



Levelling pad 

Cl 

$10/m2 , 

Wall Fill 

C2 

^Z/mGKg 

Cost of Geotextile 


$(7;(0.03) + 2.6)/m2 

Cost of Geogrid 


$(T^(0.04) + 2 Q)/m^ 

Cost of MCU face units 

C4 

%0/rrF 

Engineering and Testing Cost: 



Geotextilo wall 


$30/m^ 

Geogrid Wall 

C5jj 

SlO/m^ 

Installation Cost 

ce 

1 $50/m^ I 
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^(i “ ^ Hr 0.45 


(3 3) 

Since the wall is to be constructed with constant total length of the geosynthetic, 
accoxding to Rankine s theory, the most critical condition for pull-out occurs at the 
top where the embedment length, 1^, beyond the failure plane is a minimum. So, 
according to Fig , the minimum length of the reinforcement, 

dmirt = ifimm “h " Smin) tin 9 (3 4) 

Initial values to the design vector are assigned For the length of the reinforcement, 
the lecommended minimum value of i =: 0 is a rational initial value An initial 
value of the strength of the reinforcement is also chosen. 

Checking External Stability: 

Overturning. 

Fig. 3.3 illustrates the forces involved in the analysis The safety factor for 
overturning is evaluated by considering the moment about point A, the toe of the 
wall Since the depth of the soil in front of the wall is small, the stabvlmng passive 
piessure, is ignored. 


FS 


over 


3(7/^d + 


(3 5) 


Sliding: 

The forces involved in the sliding failure are shown in Fig. 3.4 The sliding 
resistance is calculated as. 


R^ = {'yjhd + q,)lm5 (3 6) 

where 5 is the interface friction angle between the soil and the fabric Bonoparte 
et al(19S7) and Je\vell(1990) have proposed that 5 = (2/3)(i. The force attempting 
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to cause sliding is the active thiust which can be estimated from the following 


equation, 


Pa = 




(3 7 ) 


Therefore, Factor of safety for base sliding is the ratio of the resisting force to 
the force causing sliding. 


FSm = ^ (3 8 ) 

a 

Bearing Capacity Analysis. 

The ovei turning moment from the lateral pressure in the backfill causes a transfer 
of bearing pressure from the heel of the wall to the toe. Hence, the base reaction no 
longer acts along the centre line but is eccentric by a distance e.(Fig 3.3) Considering 
a unit length of the wall and taking moments about the toe, 

^ + 3g,) 

Q{vh<i+Qjl) 

As with other gravity retaining walls, the middle third rule must be complied 
with in order to avoid any tension between the wall base and foundation soil Hence, 
the trapezoidal bearing pressure distribution gave the maximum bearing pressure, 

9m«»i ^ 


qmax - {jfhd + (3 10 ) 

The net ultimate bearing capacity was assessed using Terzaghi’s equation for a 
strip footing, 


Qnet - 

The safety factor for bearing capacity is given by, 
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F Shear = 


(3 12) 


Internal Stability: 


Safe design strength of the geosynthetic is equated by 


T 

T ^ 


The lateral load to be resisted by the reinforcement 


^lat !^d T 


(3.13) 


Assuming 100% geosynthetic coverage in the plan view, the geosynthetic should 
safely cairy this load. 

The approximate number of layers of the geosynthetic required is 


The number obtained is rounded and an additional layer is added to account for 
practical layout consideration. Assuming uniform spacing between the layers, the 
spacing was calculated as, 


h(l he 


Chock for Pullout: 


The factor of safety against pullout is given by the equation 


2(7;z + ?i) tan^ir^ 


( 3 . 14 ) 


(3 15 ) 


where, 

z = depth of the layer being designed. 

Pr = angle of frictional sliding resistance between the soil and the reinforcement. 
le = embedment length of each layer of reinforcement beyond the failure plane 
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TJ IS tlie tensile strength required at each Ie\el to resist the internal lateral 
picssiircs 


Ti = SfX (7/i 


wheie, 

5, — l/2(distance to reinforcing layer above + distance to reinforcing layer 
below). 

c/i = horizontal earth pressure at the middle of the layer 


< 7 * = KaXQs+lfbZ) 

Since the design Is assumed for equal spacing between the reinforcements, 

z = 2 X s 


where i vaiies from 1 to n(. 

The embedment length of each layer of reinforcement beyond the failure plane 
IS calculated as, 

I — [hd - z) tan d 

In addition, for geotextile wrap-wall , a minimum overlap length is required for 
re-embeclmcnt. the total length of the geotextile is then calculated as, 

hoi = ^ + s -h lolmtn 

In addition to static forces, the reinforced earth retaining wall may be subjected to 
dynamic forces, the following section briefly desenbes the procedures to be adopted 
in the design in such conditions. 

I Seismic Stability 

During an earthquake, a reinforced soil wall is subjected to a dynamic soil thrust at 
the back of the reinforced zone and to inertial forces within the reinforced zone in 
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addition to the normal static forces Current design procedures use a pseudo-static 
approach since the dynamic response of even the simplest type of retaining walls 
IS quite complex The common approach involves estimating the loads imposed on 
the wall during earthquake shaking and then ensunng that the wall can resist these 
loads An attempt was made to find the optimum cost of the reinforced wall when 
subjected to earthquake loading. The design steps are the same as for static analysis 
except that the earthquake loading is represented pseudostatically by the dynamic 
soil thiust, and the inertial force on the reinforced zone, Pm (Fig. 3 5). The 
external stability of the wall can be analysed by the following procedure. 

• The peak accleration coefficient at the centroid of the reinforced zone is cal- 
culated from the given value of peak horizontal ground surface accleration 
coefTicient, cr,nai, 

Olc (1.45 0^mar)0!rnos 

• The dynamic soil thrust is calculated from 

SPag = O.Z75ae7bh] 

• The inertial foice acting on the reinforced zone is 

Pm = oijfhdl 

• 6Pa, and 50% of Pm is added to the static forces acting on the reinforced zone 
and the usual checks for external stability are done The reduced value of Pm 
is taken to allow for the fact that the maximum values of Pas and Pm are 
unlikely to occur at the same time. 

The internal stability is calculated using the following steps- 

« The pseudostatic inertial force acting on the potentially unstable failure zone 
is determined by, 

Pi A ~ 
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wher Wa is the weight of the failure mass (mass contained within the Rankine 
Failure envelope). 

• PsA is distributed to each reinforcement layer in proportion to its resistant 
area (the area of the reinforcement that extends beyond the potential internal 
failure surface). This process produces a dynamic component of tensile force 
for each layer of reinforcement. 

• The dynamic component of the tensile force is added to the static component 
to get the total tensile force for each layer of the reinforcement 

• The check for the allowable tensile strength of the reinforcement is earned out 

• Check foi pullout failure is earned out, 

3.2.3 Objective Function 

The objective function is the cost estimate for the proposed design of the geosyn- 
thetic reinforced retaining wall. The costs involved per metre length of the wall are 
as follows'. 

1. Cost of Levelling Pad =ci 

2. Cost of the wall fill = 0 ^ x 'Yf/g x hd x I 

3, Cost of the geosynthetic used=C 3 xm xi 

4, MCU^ face units cost=C 4 x /ij 

5, Engineering and Testing Cost=C 5 x hd 

6. Installation Cost=c6 x hd 
^Only for Geogrid Wall 
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The objective function is the sum of all the above costs and, thus, is a function 
of the length and strength of the reinforcement which are the design variables. 

D'^ = (dijda) 

The objective function is minimized with respect to the design vector 

F^f{D)=f{dud2) 

3.2.4 Design Constraints 

The vaiious design constraints that are to be placed on the choice of design variables 
and other design parameters were as follows: 

• The initial length of the reinforcement chosen for design is to be greater than 
the critical length calculated by Equation 3.4 thus, the first design constraint 
Is, 

9}i ~ ^miTi “* dj. ^ 0 

• The initial allowable strength of the geosynthetic is chosen to be greater than 
the minimum specified. 

9J2 ~ ~ ^2 $ 0 

• The allowable strength of the geosynthetic is not to exceed the maximum 
available. 

9)Z — ^2 ~ £ 0 

• The factor of safety for overturning calculated by the relation (Equation 3 5) 
is to be greater than the design factor of safety 

— FSoVd ” FSover < 0 


91 



The fcictor of safety for sliding calculated from the relation (Equation 3 8) is 
to be greater than the corresponding design factor of safety 

• The factor of safety for bearing capacity calculated from Equation 3 12 is to 
be greater than the design safety factor. 

FShc^ — F iSjeor ^ 0 

• The spacing between between the reinforcements obtained from the expression 
(Equation 3 14) is to be greater than the recommended minimum spacing. 


9}j ®m«n S ^ 0 


• At each level of reinforcement, the factor of safety against pullout calculated 
from Equation 3.15 is to be greater than the design safety factor 

P;, = FSpij - FSpuii < 0 

The problem could be stated as an optimization problem as follows: Find the 
decision vector Djn » such that 

F == f[Djn) 

is the minimum of E(D) subject to 

The Sequential Unconstrained Minimization Technique (Fox, 1971) was used in 
whicli interior penalty function methods was employed in conjunction with. Powell s 
multidimensional search and quadratic interpolation for minimising steps. The 
interior penalty function approach needs a feasible starting point for initiating the 
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solulion. Siicli a starting point was easy to get. Results are obtained thus for a 
IcdSiblc initial design vectoi using interior penalty function approach. A detailed 
clcseiiptioii of the optimization procedures can be found in any standard text book 
on the optimization theory. 


3.3 Results and Discussions 

3.3.1 The Developed Program 

riie /low duut for the developed optimization procedure is presented in Fig 3 6 A 
hiiel dosciiption of the subroutines used in the procedure is presented here. 

• INTPEN : This subroutine described in Fig 3.7 is the Interior Penalty Func- 
tion Method in which a new function, is constructed by augmenting a 
penalty term to the Objective function. 

^(D,n)=F(D)-r,Z~7^ 

j=i 9] \F)) 

• POWEL : This subioutme utilizes Powell’s Conjugate Direction method for 
pattoin move in order to find the minima. 

• QJ'’rT : Tins subioutine approximates the given function by a quadratic func- 
tion foi which finite minimum can exist. 

• FUN ■ The subroutine calculates the value of the objective function. 

• CONSTR • This subroutine calculates the value of the design constraints. 

• Convergence is assumed to be achieved when the chamge in the value of the 
objective function between two consecutive cycles is less than the desired accu- 
racy. The same convergence criterion is used al through wherever convergence 

is to bo chocked. 
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'1 lie e/iocUvcjicsR of tlie developed methodology for the optimum design has been 
tleiHOiistinted with an example problem. Following the standard procedure in prac- 
tice, as enunciated eailier, the cost of construction of the reinforced wall is estimated 
With the design parameters of the example problem (FHWA Manual) as the starting 
point design vector, optimal cost for the wall is obtained 

3.3.2 The Example Problem 

Wall Description: The wall is about 200m long and exposed wall height is 5m. 
Smcliargo and seismic effect are ignored. The reinforced fill is imported to the site 
at a I'.ost of $3 per lOOOKg. An effective angle of internal friction, <j), of 34° and unit 
woigiit of 20KN/m^ is assumed for the fill. The retained backfill is assumed to have 
(f) = 30° and 7 <, = 19° 

Factors of Safety 

The design factors of safety chosen are, 

For external stability, 

• Sliding = 1.5 

• Overtuining = 2 

• lloariug Capacity ^ 2.0 

For Intel nal Stability, 

Reinfoicemcnt failure = 1.5. 

Pullout failure = 1.5. 

Calculations: 

The Active earth pressure Coefficients are calculated from Equations 3 1 and 3.2. 
The values are 

Ka, = 0.283 
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Ka, = 0.333 

Design height, li^ = 5 45m, from Eq. 3.3 
Assume l/h ratio of 0 7 so, Z = 0.7 x 5.45 

I = 3 815 » 4m 

Cheeks for external stability using equations 3.5,3 8 and 3.12 gave the following 
values of factois of safety. 

FSo,er = 5 11 
FSsi^d = 1 91 
FSbear = 8 62 

All those values are above the design factor of safety values, so the wall fulfills 
exteinal stability requirements. 

The Lateral Load to be resisted by the geogrid is: 

= i(0 283)(20)(5.45)2 = 84 06KN/m 

Cl Ci 

Assuming that a georid with a long term design strength of 2^KNjm is chosen, the 
sale design strength of the geogiid is given by : 

Td = Ta/1.5 = 13.33 

Tlie apjiroximato number of layers of the geogrid required are 

?«®= 6.2 

13 33 

This iiiunber is approximated to 8 to account for practical layout considerations and 
5.1 slope. Spacing of the layers is calculated by Eq 3.14 

sp - 5/8 = 0,625m 

COST ESTIMATE: 
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1 Cost of the Levelling Pad= (200m)($l0/m) = $2000 

2 Cost of the reinfoiced wall fill. 

(2007«)(5 45m)(4m)(20/irAf/m^)/(9.8) x $ 3 / 1000/^5 = $26670 

3 Cost of Gcogiid soil Reinforcement* 

8 layers(4m)(200m)($5/m^) = $32000 

4. Cost of MCU face units: 

(200m) (5 45m)($60/m2) = $65400 

5. Engineering and Testing Cost: 

(200m) (5.45m) (SlO/m^) = $10900 


6. Installation Cost: 


(200m)(5.45m)($50/m^) = $54500 


The total cost estiniate=$191470 

Cost estimate from Optimization Program == $176536 

Total Net Saving = $14934 
The percentage saving=7.79%. 
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Figuie 3 4* Forces for External Stability Analysis (Sliding) 



(n) 




Cb) 


FifiUic 3.5: Seismic Stability Analysis •(a)Geometry and notation for reinforced soil walls, (b) 
Static and pscudostatic forces acting on reinforced zone 
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Rcsfici Inipuit. 



Figure 3.G: Flow Chart for the Developed Program 


I 

i 


100 


























3.3.3 Optimum Cost Tables 


Typical optiiiiuiii cost tables for different wall heights(3 to 10m) are prepared for 
design paiainetcis as specified in Table 3 1 The walls are designed for different 
combinations of static and dynamic loading conditions. The results are presented in 
Tables 3.2 to 3.9 for both geotextile Wrap wall and Geognd wall The tables show 
that the optimum costs increase for increasing wall heights The optimum design 
vectoi for a particular combination of loadings and wall height is almost the same 
for both the types of walls sirice they have essentially the same design procedure 
The geogiid type wall shows more costs because of the additional cost of the MCU 
face units. Foi example, for a 7m high wall, the percent increase in cost for a loading 
of IQKN'/m is 10.2% over the no-loading case, and for a loading of 2QKNfmy the 
conespouding peicent increase is around 21%. The general observation is that for 
the same height of the wall and same loading condition, the geognd wall shows 
higher costs in the range of 20 — 30% for low wall heights (3m to 6m), while the 
increase is 10 - 20% for heights greater than 6m, 



lablc 3 2 Optimum Cost Table for Geotextile-wrap Wall, = 0, 5, = 0 


Ht. 

(m) 

1 

(m) 

(JCiV/m) 

Fill 

Vol 

geotextile 
area (m^/m) 

Tlf 

spacing 

(m) 

Cost 

$/m^ 

3 

2.69 

48 38 

9.29 

10:39 

2 

1.5 

77.73 

4 

3 21 

40 25 

14.3 

16 64 

3 

1 33 

92 87 

5 

3 73 

40 24 

20.35 

23 94 

4 

125 

108 37 

6 

4.25 

56.37 

27.44 

27 02 

4 

1 5 

112 49 

7 

4.78 

45,12 

35.57 

41.65 

6 

1.67 

139 67 

8 

5.29 

41.46 

44 74 

58 36 

8 

1.0 

167.1 

9 

5.84 

45 38 

55 21 

70.59 

9 

1 0 

183 19 

10 

6.49 

56.02 

68 16 

77 43 

9 

1 11 

188 75 


Table 3.3; Optimum Cost Table for Geotextile-wrap Wall, ah = 0, 17, = \GKNjm 



1 

(m) 



geotextile 
area (7n^/m) 

m 

spacing 

(m) 

Cost 

$/m^ 

3 

2.69 

31.20 

9.29 

14 09 

3 

IQ 

90.04 

4 

3.21 

32,86 

14.3 

20.85 

4 

1.0 

105.42 

5 

3.73 

35.72 

20.35 

28 67 

5 

1 0 

120.91 

6 

4.25 

39,06 

27.44 

37 52 

6 

10 1 

136 49 

7 

4.96 

42.65 

36.95 

48.73 

7 

1 0 

153 94 

8 

5.29 

40.57 

44.74 

64 65 ' 

9 ' 

0 89 

179 64 

9 

5.82 

44.60 

54.96 

7715 


09 

195.48 

10 

6.49 

44.82 

68.16 

92.41 

11 

0 91 

21197 


103 








T.ible 3.4. Optimum Cost Table for Geotext.ie-wrap Wall, a/. = 0, 9, = 2^KNIm 


Ht 

(m) 

1 

(m) 

{KN^m) 

Fill 

Vol 

geotextile 

area 

ni 

spacing 

(m) 

Cost 

3 

2 72 

r^.oo 

9 39 

17 88 


i i 

0.75 

103 79 

4 

3 37 

30.00 

14.99 

25.85 

5 

08 

119 78 

5 

4 01 

41.26 

21 84 

30 04 

5 

1.0 

124 23 

G 

4.64 

44 31 

29 94 

39.85 

6 

1 0 

140 95 

7 

5 27 

40.88 

39 28 

5718 

8 

0.88 

169.66 

8 

5.90 

44.87 

49 85 

70.10 

9 

0 89 

186.54 

9 

6.53 

48 87 

61 67 1 

84 26 

10 

0.9 

203 54 

10 

718 

48.50 

75 40 ' 

108 18 

12 

0.83 

232 72 


Table 3 5. Optimum Cost Table for Geotextile-wrap Wall, a* = 0 05, 9^=0 


Ht. 

(uv) 



Fill 

Vol. 

geotextile 
area (m^/m) 

nj 

Spacing 

(m) 

Cost 

S/m^ 

3 

2.88 

30.00 

"9.96 

14,66 

3 

1.00 

91.64 

4 

3.72 

45.66 

16.56 

18,16 

3 

1.33 

97.73 

5 

4.55 

45.66 

24.84 

27,23 

4 

1 25 

BMI 

6 

5.39 

47.96 

34.79 

37 97 

5 

1 2 

134 93 

7 

6.23 

51.19 

46.41 


6 

1 17 

153.86 

8 

7.07 

54.88 

59.72 

64,46 

7 

1 14 


9 

7.90 

58.83 

74.68 


8 

1 125 


10 

8.78 

56.49 

92.20 


10 
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Table 3.6 Optimum Cost Table for Geogrid Wall, Qh = 0, q, = 0 


Ht. 

(m) 

1 

(in) 

Tci 

(KN/m) 

Fill 

Vol 

geotextile 

area 

Ui 

spacing 1 
(m) 

Cost 

S/m^ 

3 

2 69 

48.38 

9 29 

5.39 

2 

1.5 


4 

3.21 


14 3 

9,64 

3 

1 33 

119 59 

5 

3.73 

40.25 

20 35 

14.94 

4 

1 25 

134 45 

6 

4.25 

56.37 

27.44 

17 02 

4 

1 5 

138 73 

7 

4.77 

45 12 

35.57 

28 64 

6 

1.67 

1 

165.12 

8 

5.29 

41.46 

44 74 

42.36 

1 

8 

1 0 

191 89 

i 9 

5 84 

45 38 

55.21 

52 58 

9 

10 


mm 

6 49 


68.17 

1 

58 43 

9 

1 11 

214 39 


Table 3 7. Optimum Cost Table for Geogrid Wall, = 0, q, = lOKNfm 


l-U. 

(m) 

1 

(m) 

^0 

(KiV/m) 

■ Fill 
Vol. 

geotextile 

area 

m 

spacing 

(m) 

Cost 

S/m^ 

3 

2.69 

31.20 

9.29 

8.08 

3 

1,0 

116,87 

4 

3,21 

32.86 

14.3 

12 85 

4 

10 

131 37 

5 

3.73 

35.72 

20.35 

18.67 

5 

10 

146.33 

6 

4.25 

39.07 

27.44 

25.52 

6 

10 

161.61 

7 

4.96 

42.65 

36 95 

34 72 

7 

10 

178,94 

8 

5.58 

40.57 

47.16 

50 22 

9 

0.89 

206.89 

9 

5.82 

44.60 

54.96 

58.15 

10 

09 

220.02 

1 10 

6.49 

44.82 

6816 

71 41 

11 

0 91 

J. 

236.56 
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3.3.4 Conclusions 


Optimal cofjt tables aie presented for different combinations of loading The Geognd 
walls show greater optimal costs because of the inherent design involving MCU face 
units The saving is of the order of 7 — 8% for the optimum design procedure 
developed here over the standard design 
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Chapter 4 

Scope of Future Studies 

The piosent study, with regards to reinforced foundation beds, may be further 
extended to consider- 

• The eflect of non-linearity of both soil and geotextile on the settlement and 
stress response. 

• The bending stiffness of the reinforcement m the analysis. 

• Furthei icfining of the mesh to study the effect on the stresses 

• Higher order elements to increase the computational efficiency 

As for the reinforced earth walls, future works may involve one of the following. 

• Optimal design of other reinforcing structures. 

• Finite Element Analysis of Reinforced Earth Walls 
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